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Description 

[00011 La presente invention concerne d'une manlere g6n6rale te trartement tf Image numerique pour ameliorer I'as- 
pect visuel de d'images et ainsi en faciliter {'exploitation ulterieure. 

[00021 L'Inventlon a trait plus particuliferement k un procddd de correction spatiale d'irrtages numdriques capt6es par 
des dispositifs embarqu6s dans des engins mobiles tels qu'avions ou navires. Un tel disposltif -photograph ie" une 
image ligne par ligne et reconstrtue I'image par la juxtaposition des llgnes. Or I'engin mobile subit au cours de son 
d§placement des mouvements perturbants tels que roulis ou derive (lacet) qui provoquent des distorsions de I'image 
[00031 Si le disposltit embarqu^ est fixe par rapport h Tengin mobile, un gyrom6tre de mesure de Vitesse angulaire 
peut §tre associ6 au dispositif de manlere k mesurer en temps rdel I'angle des mouvements perturbants par rapport 
k la trajectoire de I'engin et k corriger en consequence la distorsion de I'image. 

[00041 Altemativement. le dispositif embarqu6 peut dtre plac6 sur une plate-forme stabilis6e qui conige tout mou- 
vement perturbant. Le dispositif embarque ne subit alors plus de mouvement perturbant. 

[00051 Ces mat6riels de correction, tels que gyrom^tre de mesure de Vitesse angulaire ou plate-fomie stabilis6e 
sont trds coOteux et ne corrigent pas des distorsions supplementaires, sous forme d'instabilitd dans le sens ligne de 
I'image. provoqu6es par du bruit sur un signal de synchronisation Cgigue") lors de la transmission des images par 
exemple de I'engjn mobile, tel qu'avion ou navire, vers une station terrestre de trartement et diffusion ou reproductbn 
d'image. ^ 

[00061 La demande de brevet GB-A-2 168 874 divulgue un proc6d6 de correction spatiale d'une Image num6rique 
dans laqueile les lignes de pixel sont decal6es les unes par rapport aux autres principalement k cause du roulis de 
I'engin, tel qu'avion. La correction consiste k comparer deux k deux des lignes d-image perpendiculaires k la direction 
de i-avion, en accumulant les differences entre des pixels de deux lignes pour cinq d^calages de pixel et en determinant 
la difference accumuiee minimale afin rfen d6duire le decalage entre les deux lignes. Toutefois, ce decalage exprime 
en nombre entier de pixels ne refiete pas exactement le decalage reel des lignes d'Image. 

[00071 La presente invention vise k foumir un proc6d6 de traitement spatial d'Image numerique qui ameiiore I'aspect 
visuel d'une image, en corrigeant plus precisement non seulement des distorsions d'Image dues k des mouvements 
perturbants, mais egalemeht des distorsions dues k des defauts de synchronisation. 

[0008] A cette fin. un procede de correction spatiale d'image numerique composee d'un ensemble de rangees de 
pixels qui sont des llgnes ou des colonnes de IMmage. 

comprenant la comparaison des rang6es deux k deux pour determiner des premiers decalages entre rangees 
exprimes en nombre entier de pixel, 
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est caracterise en ce qu'il comprend 

la determinatbn de seconds decalages entre rangees exprimes en nombre non entier de pixel, en fonctlon des 
premiers decalages, 

la determination pour chaque rangee d'un troisieme decalage de ladite rangee par rapport k une rangee de refe- 
rence en fonction des seconds decalages, et 

la correction de chaque rangee par un decalage des pixels de ladite chaque rangee par rapport k la rangee de 
reference egal k la partie entiere du troisieme decalage pour ladite chaque rangee, et par un lissage spatial de 
pixels de ladite chaque rangee avec des coefficients de lissage fonctlon du troisieme decalage pour ladite chaque 
rangee. 



[0009] La comparaison des rangees deux a deux comprend pour deux rangees donnees. le calcul de sommes de 
correlation entre les deux rangees decaiees respectivement lorsque I'une est soumlse par rapport k I'autre k des 
decalages predetermines exprimes en nombre entier de pixel, la recherche d'une somme de correiatbn minimale 
pamni les sommes de correlation caicuiees afin d-en deduire ledit premier decalage entre les deux rangees exprime 
en nombre entier de pixel. 

so [0010] Une somme de correlation entre deux rangees decaiees d'un decalage exprime en nombre entier de pixel 
est egale a la somme des valeurs absolues des differences des valeurs des pixels d'une rangee et des valeurs des 
pixels decaies dudit decalage de Pautre rangee. 

[0011] De preference, la determination du second decalage entre deux rangees exprimee en nombre non entier de 
pixel est affinee par une interpolation des sommes de correlation fonction des decalages predetermines entre les deux 
rangees. par une fonctlon polynomiale d'ordre superieur ou 6gal k 2 autour du premier decalage. le second decalage 
etant I'abscisse d'un minimum de la fonctlon d'interpolation. Par exemple. la fonction polyndmiale peut-etre une fonction 
parabolique dont le sommet constitue ledit minimum. 

[0012] La determination pour une rangee donnee d'un troisieme decalage de ladite rangee par rapport k une rangee 



2 



EP 0 629 082 B1 



de reference est rdalis^e par la somme des seconds d^calages des rangdes comprises entre la rang§e de r^f^rence 
et la rangee donn6e. 

[001 3] Salon une autre caract§ristique, les coefficients de llssage spatial d'une rang6e sont tels que leur somme est 
6gale ^ 1. Par exemple, les coefficients de lissage spatial d'une rang6e sont tels qu'un premier coefficient est 6gal d 
s . la partie ddcimale du troisl^me d^calage de la rangee et qu'un second coefficient est 6gal au compl^mentaire k^ 6u 
premier coefficient. 

[0014] De pr6f6rence. les seconds d6calages entre rangees sont remplac6s par des d6calages corrig6s fonctions 
des seconds d^calages et tels que leur valeur moyenne est nulle. Cette 6tape est fondde sur la distinction entre les 
d§calages ^ corriger dus ^ des perturbations telles que roulis et des "d^calages naturels' dans I'lmage, tels qu'une 

10 route ou un cours d'eau non parall^le ^ la trajectoire de I'engin. 

[0015] En variante, lorsque les d^calages h corriger sont dus k de la derive, les seconds ddcalages entre rangees 
sont remplac6s par des d6calages corrig6s respect ivement 6gaux k la moyenne des seconds d6calages. 
[0016] D'autres caractdristiques et avantages de la pr^sente invention apparaitront plus clairement k la lecture de 
la description suivante de plusleurs realisations pr^fdr^es de Tinventron, en reference aux dessins annexes corres- 

is pondants dans lesquels : 

la figure 1 est une vue schdmatique en perspective d'un adronef 6quip6 d'un dispositif de prise de vue; 
la figure 2 est une representation sch6matique d'un premier type de deformation d'image; 
la figure 3 est une representation schenriatique d*un second type de deformation d'image; 
la figure 4 est un bloc-diagramme des etapes princlpales d'un procede de correction d'image selon une premiere 
realisation; 

la figure 5 est un algorithme de determination de premier et second decatages entre lignes selon une premiere 
realisatbn; 

la figure 6 est un algorithme de correction de second decalage entre lignes selon la premiere realisation; 
la figure 7 est un aigoritiime de determination de Timage corrigee selon la premiere realisation; 
la figure 8 est un bloc-diagramme des etapes principales d*un procede de correction d'image selon une seconde 
realisation; 

la figure 9 est un algorithme de determination de premier et second decalages entre lignes selon la seconde 

realisatbn; 

la figure 10 est un algorithme de correction de second decalage entre lignes selon la seconde realisation; et 
la figure 11 est un algorithme de determination de t*image corrigee selon la seconde realisation. 
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[0017] En reference k la figure 1 , dans un aeronef 1 , tel qu*un avion ou un drone, est fix6 un moyen de prise de vue 
2 comprenant essentiellement un moyen de detection de preference sensible aux infrarouges, un analyseur mono- 

35 ligne et un moyen de m6moire et/ou un moyen de transmission de donnees. % 

[0018] L'a6ronef 1 sun/ole une bande rectiligne de terrain 3 de maniere k photographier cette bande. La prise de 
vue est effectuee par un balayage eiectronique ligne par ligne. Une ligne L est une 'bande eiementaire" de la bande 
3, transversale k la bande 3, et la juxtaposition d'un ensemble de lignes successives forme une image d'une portion 
de bande 3. le deplacement de Taeronef le long de la bande fournissant le balayage d'image (frame). 

40 [0019] Une ligne captee par le moyen de detection subit notamment une conversion anatogique/numerique en entree 
de I'analyseur. Une ligne est un ensemble de M pixels et un groupe de N lignes est une nrratrice de N x M pixels, formant 
une image numerique k N lignes et M colonnes. Une ligne ou une colonne de I'image constitue une rangee de pixels 
dans la matrice, et par extension une rangee de I'image au sens de la presente invention. Une ligne numerique est 
memorisee par le moyen de memoire. En variante, les niveaux de luminoslte des pixels de la ligne sont emis par le 

^5 moyen de transmission vers un recepteur associe k un moyen de memorisation et de traitement de donnees. Le terme 
•luminosite" se rapporte k la luminance pour une image en noir et blanc, ou bien k I'une des composantes primalres 
rouge, verte et bleue ou k I'un des signaux de luminance et de chrominance pour une image en couleur. 
[0020] Le dispositif de prise de vue et en particulier le moyen de detection optlque sont fixes par rapport k I'aeronef. 
Pour obtenir des images reproduisant le plus fidelement possible la scene effective de la bande 3, I'aeronef vole de 

so preference en ligne droite et k altitude constante. Cependant, en raison des conditions atmospheriques, I'aeronef subit 
en cours de vol des mouvements perturbants notamment de roulis, c'est-^-dire des petites rotations altemdes R autour 
de I'axe longitudinal de I'aeronef. et de derive, c'est-^-dire une rotation D autour d'un axe vertical, ou de lacet, c'est- 
^-dire des petites rotations alternees autour d'un axe vertical. 

[0021] Dans la figure 2 est montree schematiquement une image perturbee uniquement par le roulis de I'aeronef. 
^ La forme geometrique reelle de la scene sur la portion de bande de terrain photographiee est un rectangle fonce, sur 
un 1ond clair. L'image foumie par I'analyseur k bord de i'aeronef soumis au roulis presente une forme geometrique 
deformee. 

[0022] L'image est composee de N lignes L^ k L^ et de M colonnes k 0^^, dont une ligne et une colonne 
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sont repr6sent6es Mitre tfexemple. avecl S n < N et 1 < m s M. Le roulis provoque une d§fonnation de I'image dans 
le sens ligne, c esl-^-dire les lignes sonl d6caldes les unes par rapport aux autres. 
[0023] En pratique, une image comprend par example 1 024 lignes et 1 024 colonnes 

[0024] En reference a la figure 3. une image sch6matique est perturb^e uniquement par des derives de Pa6ronef 
pnncipalement suivant un mSme sens par rapport k to tiajectoire de ra§ronef. La forme g6om6trique rdelle est k 
noowauunrBclanglefoncdsurfondclalrLlmagefoumieparl'analyseuraborddel'aero^^^ 
rZ^f,^^^'"""* d6fomf,6e rfeuNant de la transformation du rectangle d'origine sensiblement en losange 
VT^^x L„ PanTi' N lignes et une colonne C„ parmi M colonnes sont dgalement representees dans la figure 

3. La derive provoque une deformation dans le sens colonne. c'est-S-dire les colonnes sont decaiees les unes par 
rapport aux autres. 

P02q Enpratiquejesdefoimatjonsduesauxroulisetlesdeformationsduesiladerivesecomblnentaieatoirement 
^ .TJl""^^ P0rte"ses en g6n6ral d'un grand nombre de details, ces deux types de detonration 

rendent difficiles la lecture et I'interpretation des images par un operateur. De plus, d-autres defonnalions alierent une 
.mage. En parliculier. dans le cas de la transmission tfimages entre Paeronet et un mcyen de reception, la synchroni- 
^.ondesdc^inees est realiseeparune reference deseuilsurunfrontmontant d'un sign^ 

le front rnontant est bmrte. la position du seuil sur le front montant subit des variations provoquant un phenomfene de 
gigue (jitter), c est-a<lire d'instabilite dans le sens des lignes de llmage. Cette instabilite piovoque des decalaqes entre 
lignes analogues aux deformations dues au roulis. mais de plus faibles amplitude 

^|En reference^ la figure 4. un procede de correction d'image destine plus particulierement k une image 
d6fonT)«e par du roulis est illustr6 sous la fomie d'un algorithme compose de cinq etapes prhcipales E10 k E50 et 
programme dans une memoire de calculateur. La correction d-lmage est realises dans I'aeronef. ou en variante au sol 
Dans tous les cas. comme il a d6ja ete precise, une image est une matrice de N x pixels memorls6e Le precede 
comprend une prerpiftre etape E10 de lecture d-une image qui consists a lire les niveaux de luminosite des N x M 
pixels de I image k N lignes et M colonnes dans une memoire d'image. 

[(KttBI L-etape E20 detemr<!ne et memorise les d6calages des lignes prises deux k deux successivement. L'etape 
f^f compos6e d un ensemble d-etapes E21 k E27 qui seront ddcrKes en detail en reference k la figure 5 
[0M9] L etape E30 traite les decalages entre lignes pour separer les decalagesdus au roulis. de tout decalage non 
le ^ la prise de vue. mais existant reellement sur la portion de terrain photographiee. par exemple une route formant 
un angle avec les bortfe de I'image. ou un cours d-eau qui serpente. L'etape E30 est composes d'un ensemble d-etapes 
E31 k E34 qui seront decntes en detail en reference k la figure 6. 

[OMOJ L-etape E40 detemiine rimage corrig6e k partir des decalages entre lignes d6temiln6s k l'etape E30 En 
varian e. I etape E40 detemiine I'image corrigee k partir des decalages entre lignes determines k retape E2a cette 
variante ne comprend pas l'etape E30. 

3S ^^"l f ? '^^ T- ^'^^ "^^^ msmtote d'etapes E41 k E46 qui seiont decrites en detail 

en rdTerence a la figure 7. ^ 

^"''^J!® ^"^'^ memorises k l'etape E50. et est. de preference, egalement 

afficfiee k I ecran d-une console de visualisation associee au calculateur. 

[0033] En variante. I'image corrig§e k l'etape E40 peut Stre traitee k nouveau k partir de l'etape El 0. Piusieurs cycles 
de correction successifs ameiiorent plus I'image qu'un cycle unique. Le nombre de cycles de correction est soit fixe 
r r ^' "''^'"'"^ 69al a 1 . 2 ou 3. soit choisi par Pop6rateur. soit encore determine dynami^ 

quementparJalgorithmesur un criterode qualite defirnage corrigee" , „ . ? 

K?t^^*7I'"^"'"*^***'^*'^"~*' '«<'«»*™"na«on decalages entre lignes comprend un ensemble d'etapes 
I ■Ati. V ^"f ^ ''i'"«9« comparee k la ligne adjacente I,., pr6cedemment m6moris6e. 
L etape E21 initialise la detemiination k la ligne qui est compar6e k la ligne L,. Dans la suite, la ligne L, n'est pas 
atfectee par a correction d'image et est une ligne de reference. En variante. i'orelre des calculs est inverse, la ligne de 
reference est la l.gne L^, et comparaison commence alors k la Hgne L^,. qui est c6mpar6e k la ligne L^. Seton une 
autre variante. la ligne dereference est une ligne quelconque del-image. a hm "ne 

. ^'^^^ decalage en nombre entier de pixel Dec k une vaieur -Dmin. Le decalage Dec est 

le decalage de comparaison de la ligne L„ avec la ligne L^,. 

[003q Pour une vaieur de decalage Dec donn6e. l'etape E23 compare la ligne I, & la ligne L,,. en calculant une 
somme de correlation Sp^ entre les deux lignes. pour le decalage Dec. U somme Sq^c est egare d : 
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ou : 

P'Xn-i,m ®st le niveau de luminosity du pixel au croisement de la ligne L^.i et de la colonne C^, 

,m+Dec ^® niveau de luminosite du pixel au croisement de la ligne L„ et de la colonne C^Dg^, 

Bord est un entier, strlctement infdrieur M, 6gal S un nombre de pixels situ6s sur chacun des deux bords de 
{'image, et non consid^r^s dans la comparaison; k titre d'exemple, pour des lignes de 1024 pixels, les bords non 
considdrds sont d'une largeur comprise entre 50 et 100 pixels. 

10 

[0037] La valeur de S^^ est m^moris^e, puis le ddcalage de comparaison Dec est incr6ment6 de 1 . 
[0038] A l'6lape E24, le d6calage Dec est compar6 ^ une valeur pr6d6termin6e Dmax qui est la valeur maximale du 
d^calage Dec. Si la valeur du d^calage Dec est infdrieure ou dgale ^ la valeur Dmax. I'etape E23 est k nouveau 
d6roul6e. Sinon, toutes les sommes de correlation S^^ pour Dec variant de -Dmin ^ Dmax sont d6termln6es, c'est- 

15 ^-dire des sommes au nombre de (Dmin + Dmax + 1 ). et I'algorithme passe k Tdtape E25. 

[0039] Les bomes de variation -Dmin et Dmax du d^catage Dec sont k titre d'exemple -25 et +25 pixels. 
[0040] L'dtape E25 recherche la somme minimale Soe^jp parmi les sommes d^tenmin^es pr^^demment. ainsl que 
le d^calage Decmin associd. Pour trouver la somme minimale Sp^^;„ parmi (Dmin + Dmax + 1) sommes, les vateurs 
adjacentes k Sq^ sont compar6es deux k deux. La somme minimale Spec^ip est la premidre somme trouv6e dont les 

20 deux sommes adjacentes pr^c^dente et suivante tui sont sup^rleures lorsque Dec croit de - Dmin k Dmax. En variante. 
toutes les sommes Sqqc sont analys^es, en comparant les sommes deux k deux, et en gardant comme r^suttat de 
chaque comparaison la somme Sq^c plus f alble. Selon une autre variante, la somme minimale Specmin premiere 
somme trouv6e dont les deux sommes adjacentes pr^idente et suivante lui sont sup^rieures de valeurs de difference 
sup6rieures k un seuil predetermine. 

2S [0041] Dans tous les cas, k la somme minimale trouv6e SQ^cmin associ6 le d6calage de la ligne L^ par rapport 
k la ligne Lp.^ exprimd en nombre entier de pixel Decmin. Or, un ddcalage exprim6 en nombre entier de pixel ne sert 
k corriger une image que d'un d^calage relatlvement 6\6v6 d'une ligne k I'autre. En praitique, le d6calage d'une ligne 
k Tautre survenu au cours de la prise de vue est fr^quemment Inf^rieur au pixel; I'tnvention pr^voit alors d'affiner le 
calcul du d^calage. 

30 [0042] L'6tape E26 determine le d6calage D^ de la ligne L„ par rapport k la ligne L^.^ , avec une precision inf6rieure 
au pixel; d'une mani6re g6n6raie, D^ est un nombre non entier de pixel. Le d6calage D^ est determine en fonction des 
ddcalages Decmin, (Decmin*1 ) et (Decmin + 1 ) ainsi que des sommes de correlation associ^es k chacun de ces 
decalages Doecmin > ^oecmin-i ^Decmin>i- ^Si^^ un rep^re ayant les d^catages pour abscisses x et les sommes de 
correlation pour ordonndes y, ces trois couples de d^calage et somme definissent trois points. En recherchant la 

35 parabole d'6quation : 

2 

y = ax +bx + c 

40 passant par ces trois points, est d6duit le sommet de cette parabole qui a pour ordonn6e la somme minimum 

affinee et pour abscisse le decalage plus prdcis D^ de la ligne L„ par rapport k la ligne L^.^. Le d6calage D^ r6sulte de 
la resolution du systeme k trois equations de parabole correspondant aux trois points afin d'en deduire les coefficients 
a, b et c. Ce decalage est ainsi egal k : 

- b/(2a) = Decmin + 0.5x{Sd^5„ ^ , - 

^Decmin-1 V X Sq^j^ - ^Decmln-1 " ® Decmin + 1 ) 

50 [0043] La valeur de D^ est ensuite m6moris6e et le paramdtre n est incremente de 1 . . 

[0044] En variante, le decalage Dp peut etre determine comme I'abscisse du minimum de la parabole passant par 
les trois.points (Decmin, So^cmin). (Decmin + k, SQec^jp^ k). (Decmin-k, So^n^j^j^) ou k est un entier. valant 2 ou 3 par 
exemple. Selon une autre variante, la parabole est remplacee par une fonction polynomialed'ordre plus eieve que deux. 
[0045] L'etape E27 compare la valeur de n avec le nombre total de lignes N. Si n est inf6rieur ou 6gal k N. il reste 

55 au moins une ligne pour laquelle le decalage D„ n'a pas ete determine et le calcul reprend pour une nouvelle ligne k 
partir de l'etape E22. Sinon, tous les decalages D2 k D^^ sont determines et I'algorithme passe k l'etape E31 qui est 
decrite dans la suite. 

[0046] En variante. une ligne quelconque L^ n'est pas comparee la ligne adjacente L„.^. nnais k la tigne i^i, ou i 
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hI!"^ ^'^ ^ ^ Para^^tre n est incr6mente de i. A tilre d'exemple. rentier i vaut 2 et seules les 

Jgnes de rang pair sonl trait^es. Les lignes de rang impair sont d6cal6es (6tape E40) d-un decalage dgal au d^lage 

rJ,Snt!fr« H f •• f^^* °" dependent des d6calages dus au roulis, mals ^galement de -dficalages' 
une oute n<« paralfele aux bords de la bande de terrain. La correction des d6calages repose sur le c^Sal que 

aSril^Sf IT"' '^y^"" ^« «9"««' 1024 p^r exemple. ta^« 

™ln.o f ^ i^'" «^ ^» ^" P'^^' o-^Jre sur les dScalages- 

en vanante, la correction est d'ordre plus 6lev6. wowaiaaco. 

f°°*? . ';f foffection des d6calages est telle qu'un ensemble pr6d6temiin6 de lignes de I'image n'est pas modifid 
lors de a dSlerrn.nat «n de I'image corrig6e. qui sera d6crite en detail dans la sufte A titre d'exemple sVSeTN 

Zlre v«rZ!T 1^°?; ^''"^'^ ^'^ °" ^' r6guli&rement r6parties ou non dans rimage. Selon 

f^i« , !'°'' '^^^'n^'^"^^ S„ S2 et P utilises pour la correction des d^calages. Le .LfRclent 

2S2yN2relatrfSla,«nat>ondudecalage des lignes par rapportabpremlfereligneder«6^^^^ 

n *!?'!f f ^ '"^^^ '^ dStemiination de la correction de d^calage k la deuxifeme ligne puisque. comma » a 

Esf its ifuf JfrT' T ^'J'*^"^" "^'"'^ determination de I'mage Srrig6e. 

E33 calcule un coefficient correcteur Cor„, qui est 6gal k S1/N + P x (n - N/2) afin de ramener la variation du dScaboe 

du dSTn ' If"""'" '''''^"'^^ "-e coefficient Cor„ esfe^ut 

'^^^ "^^^ Le paramWre n est incr6ment6 de 1 . 
H^IL J^ » '^'^ 'f f^""^^'^ " ^ valeur N; si n est inf6rieur ou 6gal a N. I'6tape E33 est a nouveau 
d6roul6e pour cette nouvelle valeur de n. Sinon. tous les decalages corrig6s Dc„ pour 2^nsN sont d^terminT e! 
I'algoritf^me passe a l'61apeE41 qui est d6taill6e dans la suite « Po ^-n.N. sont determines et 

E^^el esfefftltr l^n^^^^^^^^^^ ''"^^^ '"'^'^^ <'«^= ^ 

f^f ^ , ^'^^^ '""'^"^ correction de Pimage ^ la deuxidme ligne Lg de I'image. A l'6tape E42 est determine 

tTli^ 1 ^" ^^'^'n ^^^'"^ decalages a pour ^isn 

en^ere est no ee En. La partie decimate de Decal„ est notee k„, et le compl6mentaire S 1 de Ic,, est rr^6 L"6tew 
E43 initialise a 1 un parametre m. representatif du rang rfun pixel dans la ligne ^' ^ 

Eoe'"^ «^li!!lifl'"^!."' " f '""""^"'^ ''"'^^ P*^"'" P««' ^« dans la ligne L„ par un 

Se SrS^Tn '"''7'" ""^ "'^Se spatial de vateurs de pixel de la ligne. Le ifssage 

[t n ll^ ir ? «dj«=«"l« de rangs m + En et m + En + 1 dans le pixel corrig6 

plZnZts ^- ^T ZlS To ''''' ^' respectivemL ponder'es 

cnp^fn^l v. . ! "''"^^ ^' P^'^'"^^'« est mcremente de 1 . En variante, 

?0 5 Se on un^^ sont determines differemmenl, et sont par exemple egaux 0 et 1 respectivement. ou a 05 

„ , 7^ ^'.xHEn-i. "^.m+En P'Xn.m+En+i- Dans tous les cas. les Coefficients de ponderation sont de Dr6fe- 

renoa tels que leur somme est egale ^ 1 . de sorte qul^ ^ 

[0055] Uetape E45 verifie si le parametre m est interieur ou egal au nombre M de cotonne de Pimage et done au 

Te Jdete'i^iSe ?I "^I' " T °" ' " - dont la vaTeur co rJeTn-a 

le Lrltrrr^Hn? ^*«Pf^"4f Parcourue. Lorsque tous les pixelsde la ligne sont corriges. alors 
l^kUEntZTl ^« 1 ff^' ^ "9"^ suivante. L'etape E46 verifie si le pai^etre n esl Inferieur 
.^.tP^Lclnr f ^ ^^^""^ deroule les etapes E42 ^ E45 pour une nouvelle ligne. Lorsque 

lou es les lignes sont corrigees. Palgorithme passe a retape E50 de memorisation de I'image corrig6e CommH a 

Se'pater."^^^^^ 

[0056J En reference k la figure 8. un precede de correction d. deformation dimag. due k la d^ive est illustre 
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sous la forme d'un algorithme compost de cinq Stapes principales E100 h E500, analogues aux stapes E10 k E50 
pr6c6demment d^crites en r6f6rence k la figure 4. La prlnclpale difference est que la correction de Timage est effectu6e 
colonne par colonne, et non ligne par ligne. A L'6tape E100 sont lus N x M niveaux de luminosity de pixel m6moris6s 
d'une image k N lignes et M colonnes. 

s [0057] L'^tape E200 determine et memorise les d^alages des colonnes prises deux d deux successivement. L'^tape 
E200 est composde d'un ensemble d'6tapes E210 k E270 qui seront dicriles en detail en rifdrence k la figure 9. 
[0058] L'6tape E300 traite les d^calages entre colonnes pour sSparer les d^calages dOs k la derive de tout d^calage 
non lid k la prise de vue, mais existant r6ellement sur la bande de terrain photographi6e. L'6tape E300 est compos6e 
d'un ensemble d'dtapes E310 k E340 qui seront ddcrites en detail en r6f6rence k la figure .10. 

10 [0059] L'6tape E400 determine I'image corrig6e k partir des d6calages entre colonnes determines k retape E300. 
En variante, I'etape E400 determine rimage corrigee k partir des decalages entre colonnes determines k I'etape E200; 
cette variante ne comprend pas retape E300. 

[0060] Dans tous les cas, retape E400 se compose d'un ensemble d'etapes E410 k E460 qui seront decrites en 
detail en reference k la figure 11 . 

75 [0061] La nouvelle image corrigee des effets de la derive est memorises ^ retape E500 et est, de preference ega- 
lement affichee k I'ecran d'une console de visualisation associee au calculateur. 
[0062] En variante, rimage corrigee k i'etape E400 peut dtre traitee k nouveau k partir de retape E100. 
[0063] Comma montree k la figure 9, la determination des decalages entre cobnnes comprend les etapes E210 k 
E270, qui sont respect ivement analogues aux etapes E21 k E27 decrites en reference k la figure 5, k la difference 

^0 pfks que les decalages ne sont pas determinees ligne a ligne, mais colonne k colonne. 

[0064] La lettre 'C est ajoutee aux noms des differents parametres. Les colonnes sont comparees deux k deux 
pour determiner un premier decalage DecCmin de ta colonne par rapport k la colonne C,^.^ , exprime en nombre 
entrer de pixel, puis un second decalage DC^^ de la colonne C„ par rapport k la colonne C,^^ exprime en nombre non 
entier de pixel. Lorsque tous les decalages DC^. pour 2<m<M, sont determines, I'algorithme passe k I'etape E310. 

25 [0065] En reference k ta figure 10, le traitement des decalages entre cotonnes comprend quatre etapes E310 k 
E340. Ce traitement a pour but de corriger les valeurs des decalages DC^, pour 2^m<l^, determines k I'etape E260. 
[0066] L'etape E310 comprend le calcul de la somme SC et de la moyenne MO des decalages 00^^, pour 2^m^M. 
[0067] L'etape E320 initialise le parametre m ^ 2 et I'etape E330 determine le decalage corrige DCc„ de la colonne 
C^, qui est egal k la moyenne MC des decalages DC„. Le paramdtre m est incremente de 1 pour passer k ta colonne 

30 suivante. A retape E340, si m est interieur ou egal k M, l'etape E330 est k nouveau derouiee pour cette valeur de m; 
stnon, I'algorithme passe k retape suivante E340. 

[0068] Cette correction de derive suppose que la derive de I'aeronef est quasiment constante pendant la prise de 
vue de I'image ce qui implique que la correction est identique pour toutes les colonnes. 

[0069] En variante lorsque la derive de I'aeronef est modifies graduellement pendant la prise de vue de I'image. le 
35 calcul de la somme SO et de la moyenne MO des decalages DC^, pour 2<m<M, est remplace par le calcul d'un 
ensemble de sommes et moyennes calcuiees sur des groupes de valeurs de decalage. Par exemple, lorsque la derive 
presente trois valeurs successives pendant la prise de vue, les sommes et moyennes sont: 



40 



ml 



SCi = £ DC^ ; MCi = SCi/ml 
m 2 



45 



SO 



m2 



SC2 = X ^^m ; MC2 = SC2/(in2-ml) 

m ml- 1 



55 



SC3 = £ DCni ; MC3 = SC3/(M-m2) 

m 17121 
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avec 2<mian2<M. Le d6calage corrige DCc„ est alors dgal soft k SC,. soft a SC^ soft k SC, suivant que m est 
compreentre2etm1,oum1 + le«m2,oum2 + 1etM. suivani que m est 

[Sk ^" ■ comprend pas les 6tapes E310 k E340 pr6cMemment dfcrites. Dans ce cas 

lalgorrthmepassedirectementdel'6tapeE270&l*etapeE410. t^a.« i-e cas. 

I0071J En r6f6rence ^ la figure 11 . la determination de I'image corrig6e comprend les etapes E410 a E460 resoec- 

n est pas effectu6e hgne par ligne, mais colonne par colonne. Pour chaque cdonne C„ sont d6termin6s les niveaux 
de lum.nos.t4 corng^ de pbcel PC,„. pour 1^.N, chacun en fonction de niveaux enmgistr^s ou^SliT*. S 
l^ECm.m ei Kix^Ecm*i.ro °o deux pixels adjacents dans la mSme colonne. Lorsque toutes les colonnes sont 
corrig^es. ralgorithme passe k Tfitape E500 de memorisation de rimage corrig6e 

[0072] Les deux types de correction d-image respectivement suivant les 6tapes E10 a E40 et E100 a E400 sont oar 
deneffectuer que run des deux types de correction pour une Image donn6e. 



Revandleatlone 



1. Proc6d§ de correction spatiale d-image num6rique compos6e d'un ensemble de rang6es de pixels (L, a U" C, & 
Cm) qui sont des lignes ou des colonnes de Pimage. comprenani * 1 H<. 1 « 

r^^"" ^ (E25.E250) des premiers decalages 

entrerang6es(Decmin;DecCmin)exprim6s en nombreentierde pixel, 

caracterisd en ce qu'il comprend 

la determination (E26,E260) de seconds decalages entre rangees (D„:DC„) exprimes en nombre non entier 
de pixel, en tonction des premiers decalages. n m/ r 

SnS!T'"r'T <^^;E42j))P°"^'^f«q"« ranges (L„;C„) d'un troisiSme decalage (Decal„;DecalCJde ladite 
rangee(yc )parrapportaunerang6eder6f6rence(^,C,)enfonctiondessecondsd6.;klages(D ^ ) et 
te correcton de dTaque rang6e (E44;E440) par un decalage des pixels de ladfte chaque rangee par raSport 

hssage spatal des pixels de ladite chaque rangee avec des coefficients de lissage (k„,.k^ kC„, kC ,) fonc- 
tion du troisieme decalage pour ladfte chaque rangee. lN.1.'W.'^^m1.Kl'm2)I0nC 

* revendication 1, caracterise en ce que la comparaison (E23,E230) des ran- 
gees deux k deux comprend pour deux rang6es donnees (L,.L„.,;C„,C^,), le calcul de sommes de correction 
S^dSi^i llT. deux rangees decaiees respectivement lorsque I'une est soumise par rapport k I'autre k 

rt^i^l^n n <fo«7-T®D«»Cn*.) P^""' "^s sommes de correlation calcul6es af in d'en deduire ledit premier 
decalage (Decmin:DecCmin) entre les deux rangees exprime en nombre entier de pixel. 

3. Precede de correctkw conforMe k la revendication 2. dans lequel une somme de correlation (Sn Sn^) entre 
deux rangees decaiees d'un decabge (Dec;DecC) exprime en nombre entier de pixel est e^ale^b ^me det 

decalage (Dec;DecC) de Tautre rang6e. 

uS,f H '°''=°''''' ^''P""'^^ ^" de pixel est realisee pa 

tesdeux rangees, par une fonction polyn6mBled'ordresuperieurouegald2a^^^^ 

DecCmin). le second decalage (D„;DC„) etant I'abscisse d'un minimum de la fonction d'interpolatk.n 

;. Pr^ede de correctfon conforme k I'une quelconque des revendteations 1 a 4, dans lequel la determination (E42- 
E420) (XHir une rangee donnee (L„:C J d'un troisieme decalage (Decal„;DecalC J de ladfte rangee (L^ Vp5 

rapporlSunerangeedereference(L,:C,)estrealiseeparlasommedess^dsd2^abges(D„;Dl!)^^^ 
comprises entre la rangee de reference (L,;Ci ) et la rangee donnee (i^;C J. *" ^ 

. Precede de correction confomie k Tune quelconque des revendicatkxis 1 a 5, dans lequel les coefficients de 
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lissage spatial (k„i,k„2;kC^i,kC„2) d'une rang6e (L„;Cn,) sont tels que !eur somme est 6gale ^ 1. 

Procddd de correction conforme k Tune queiconque des revendications 1^6. dans (equal les coefficients de 
lissage spatial (k^^ K2^^^mv^^m2i ^'^^^ rang^e (L^;Cfn) sont tels qu'un premier coefficient (k„i;kCn,^) est 6gal 
k la partie decimate du troisi^me dScalage (Decal„;DecalCfyj) de la rang6e (L^;Cf„) et qu'un second coefficient (k„2* 
kC^2) ^9^' comptSmentaire ^ 1 du premier coefficient. 

Proc^d de correction conforme k Tune queiconque des revendications 1 k 7, dans lequei les seconds d^alages 
entre rangdes (D^) sont remplac^s (E33) par des d^calages corrigds (DCf^) qui dependent des seconds ddcatages 
et dont leur valeur moyenne est nulle. 

Proc6d6 de correction spatiale d'image conforme k I'une queiconque des revendications 1 k 7, dans lequei les 
seconds ddcalages entre rang^es (DC„) sont remplac^s (E330) par des ddcalages corriges respectivement 6gaux 
k la moyenne (MC) des seconds ddcalages. 



Patentanspruche 

1. Verfahren zur raumlichen Korrektur eines digitalen Bildes, bestehend aus einer Gesamtheit von Reihen von PIxeIn 
20 (L<, bis L^; bis C^). die Zeilen Oder Spalten des Bildes sind, umfassehd 

den paarweisen Vergleich (E23;E230) der Reihen, urn erste Verschiebungen zwischen Reihen (Decmin;Dec- 
Cmln). ausgedruckt in einer ganzen Zahl von Pixein, zu bestimmen (E25;E250), 

25 dadurch gekennzeichnet. dass es umfasst 

die Bestimmung (E26;E260) zweiter Verschiebungen zwischen Reihen (DniDC^), ausgedruckt in nichtganz< 
zahlrgen Pixelzahlen, in Abhangigkeit von den ersten Verschiebungen, 

fur jede Reihe (LniC^) die Bestimmung (E42;E420) einer dritten Verschiebung (DecalniDecalC^) der genann- 
30 ten Reihe (Ln;C„,) relativ zu einer Bezugsreihe {li»0^) in Abhangigkeit von den zweiten Verschiebungen (D^; 

DC J. und 

die Korrektur jeder Reihe (E44;E440) durch eine Verschiebung der Pixel dieser Re^e relativ zur Bezugsreihe 
um den ganzzahligen Anteil der drrtten Verschiebung f Or diese Reihe und durch eine raumllche Glattung der 
Pixel dieser Reihe mrt Glattkoeffizienten (k^i .kn2;kCn,i ,kC„2) Abhangigkeit von der dritten Verschiebung fur 

35 diese Reihe. 

2. Korrekturverfahren nach Patentanspruch 1 , dadurch gekennzek^hnet, dass der paarweise Vergleich (E23;E230) 
der Reihen fOrzwei gegebene Reihen (Ln.L„.^;C^,Cn^i) die Berechnung von Korrelationssummen (SoeclSoecc) 
jeweils zwischen den beiden verschobenen Reihen umfasst, wenn die eine relativ zur anderen vorbestimmten 

40 Verschiebungen, ausgedruckt in einer ganzen Pixelzahl (Dec;DecC), unterliegt, die Suche nach einer mlnlmalen 

Korrelationssumme (SD^^i^iSD^jjCmin) unterden berechneten Korrelationssummen. um darausdie genannte erste 
Verschiebung (Decmin;DecCmin) zwischen den beiden Reihen. ausgedruckt in einer ganzen Pixelzahl, abzulelten. 

3. Korrektun/erfahren nach Patentanspruch 2, in dem eine Korrelationssumme (SoeclSDecc) zwischen zwei, um eine 
^ Verschiebung (Dec;DecC), ausgedruckt in ganzen Pixelzahlen, verschobenen Reihen, der Summe der Absolut- 

werte der Drfferenzen der Werte der Pixel einer Reihe und der Werte der verschobenen Pixel der genannten 
Verschiebung (Dec;becC) der anderen Reihe gleich ist. 

4. Korrekturverfahren nach irgendeinem der Patentanspruche 2 oder 3, in dem die Bestimmung (E26;E260) der 
so zweiten Verschiebung zwischen zwei Reihen (D^iDCn), ausgedruckt in nichtganzen Pixelzahlen, durch eine Inter- 
polation der Korrelationssummen (SDeciSoecc) erfolgt, die Funktion der vorbestimmten Verschiebungen (Dec; 
DecC) zwischen den beiden Reihen sind, durch eine Poly nomialfunkt ion einer Ordnung grosser oder gleich 2 um 
die erste Verschiebung (Decmin;DecCmin), wobei die zweite Verschiebung (DnjOCJ die Abszisse eines Mini- 
mums der Interpolattonsfunktion ist. 

55 

5. Korrekturverfahren nach irgendeinem der Patentanspruche 1 bis 4, in dem fur eine gegebene Reihe {L„;C^) die 
Bestimmung (E42;E420) einer dritten Verschiebung (Decaln;DecalC^) der genannten Reihe (LniCn,) relativ zu 
einer Bezugsreihe (Li;C,) durch die Summe der zweiten Verschiebungen (DniDCn,) der Reihen erfolgt, die sich 
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zwischen der Bezugsreihe (li,C^) und der gegebenen Reihe (Ln;C J befinden. 

6. Korrekturverfahren nach irgendeinem der Patentanspruche 1 bis 5. in dem die raumlichen Glatlkoeffizienten {k^^, 
K»&kCmi.kC^ einer Reihe (Ln;Cm) derart sind, dass ihre Summe gleich 1 ist. 

7. Korrekturverfahren nach irgendeinem der PatentansprOche 1 bis 6. In dem die raumlichen Glattkoeffizienten {k„^. 
kn2:kC^i.kC^) einer Reihe (L„;C^) derart sind, dass ein erster Koeffizient (k^i.kC^i) 9«eich dem Dezimalantel'l 
der dritten Verschiebung (Decal„;DecalC J der Reihe (L„;C J ist und daB ein zweiter Koeffizient (k„2.kC„2) glefch 
der Erganzung zu 1 des ersten Koeffizienten ist 

8. Korrekturverfahren nach irgendeinem der Patentanspruche 1 bis 7. in dem die zwelten Verschiebungen zwischen 
Relhen (D„) durch korrlgierte Verschiebungen (Dc„) ersetzt werden (E33). die von den zweiten Verschiebungen 
abhangen und deren MIttelwert gleich Null ist. 



9. 



Verfahren zur raumlichen Blld korrektur nach irgendeinem der Patentanspruche 1 bis 7, in'dem die zweiten Ver- 
schiebungen zwischen Reihen (CXJ J durch korrlgierte Verschiebungen ersetzt werden (E330). die jeweils dem 
MIttelwert (MC) der zweiten Verschiebungen gleich sind. 



Claims 



4. 



A method for spatial correction of a digital image made up of a set of rows of pixels (L. to U C. to Cm) which are 
lines or columns of the image, comprising 

comparison (E23; E230) of the rows in pairs in order to detemnine (E25; E250) first inter-row shifts (Decmin- 
DecCmin) expressed as a whole pixel number. 

characterised Jn that It comprises 

determination (E26; E260) of second Inter-row shifts (D„; DC^) expressed as a non-whole pixel number, as 
a function of the first shifts, 

determination (£42; .E420) for each row (L„; C J of a third shift (Decal„; DecalC J of the sakJ row (L. C„) with 
respect to a reference row (Li; C^) as a function of the second shifts (D^; DC^). and 
correction of each row (E44; E440) by a shift of the pixels of said each row with respect to the reference row 
equal to the whole part of the third shift for said each row. and by a spatial smoothing of the pixels of the said 
each row with smoothing coefficients (k„i,k„2:kC„,.kC„2)asafunction of thethird shift fo^ 

A correction method according to Claim 1 . characterised in that the comparison (£23; E230) of the rows In pairs 
comprises for two given rows (1^. L^^; C„. C^^). cafculating correlation sums {S^; S^e^) between the two 
respectively shifted rows when one is subjected with respect to the other to predetermined shifts expressed as a 
whote^pixel number (Dec; DecC).^searchlng for a minimal correlation sum (Soecmin: SDecCmin) amongst the calcu- 
lated correlation sums In order to deduce therefrom the said first shift (Decmin; DecCmin) between the two rows 
expressed as a whole pixel number. 

A correction method according to Claim 2. wherein a correlation sum {So^^ So^cc) between two rows shifted by 
a Shift (Dec; DecC) expressed as a whole pbcel number Is equal to the sum of the absolute values of the differences 
of the values of the pixels of one row and the values of the pixels shifted by the said shift (Dec; DecC) of the other 
row. 



A correctron method according to either Claim 2 or Claim 3. wherein the determination (£26; E260) of the second 
shift between two rows (D„; DC„) expressed as a non-whole pixel number is effected by interpolation of the cor- 
rebtron sums (Sd^ Soeec) as a function of the predetermined shifts (Dec; DecC) between the two rows by a 
polynomial function of an order greater than or equal to 2 around the first shift (Decmin; DecCmin), the second 
shift (D„; DC J being the abscissa of a minimum of the interpolation function. 

A correction method according to any one of Claims 1 to 4. wherein the determination (E42 E420) for a given row 
(U; Cm) of a third shift (Decal„; DecalC J of the said row (1^; C„) with respect to a reference row (L,; C,) is 
embodied by the sum of the second shifts (D„; DC J of the rows contained between the reference row (L,- C,) 
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and the given row (L„; C^). 

A correction method according to any one of Claims 1 to 5. wherein the spatial smoothing coefficients (km, kn2> 
kC„^, kC^ of a row (L^; Cj„} are such that their sum is equal to 1. 

A correction method according to any one of Claims 1 to 6, wherein the spatial smoothing coefficients (k„-,, k^g; 
kC^i, kC„2) of a row (1^; Cj^) are such that a first coefficient (k^; kC^^) is equal to the decimal part of the third 
shift (Decalft; DecalC^) of the row (Ln; C^) and that a second coefficient (kf^; kCn2) ^^lual to the ones complement 
of the first coefficient. 

A correction method according to any one of Claims 1 to 7, wherein the second inter-row shifts (D^) are replaced 
(E33) by corrected shifts (Dc J which depend on the second shifts and the average value of which is zero. 

An image spatial correction method according to any one of Claims 1 to 7, wherein the second inter-row shifts 
(DCfn) are replaced (E330) by corrected shifts respectively equal to the mean (MC) of the second shifts. 
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FIG. 5 
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FIG. 6 
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FIG. 9 
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